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Liquids that form due to dynamics of the molecules that depend on the local 
density 
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RNA molecules in living cells form what look like liquid droplets formed by liquid/liquid phase separation. But 
unlike the molecules in conventional phase separating mixtures, RNA molecules are transported by molecular 
motors that consume energy and so are out of equilibrium. Motivated by this we consider what sort of 
simple rules for the dynamics of model rnRNA molecules lead to liquid/liquid phase separation. We find that 
dynamics that slow as the local density of molecules increases, drive the formation of liquids. We also look at 
the analogous separation of the two blocks of a block copolymer, in which the monomers of one block have 
dynamics that depend on the local density of monomers of that block. We find that this block condenses and 
separates from the monomers of the other block. This is a simple model of the out-of-equilibrium domain 
formation found in the chromatin in the nucleus of cells. 


The contents of living cells are in the liquid state. But 
these contents are not in one liquid phase, they appear 
to be in a number of coexisting liquid phases. There are 
what look like liquid droplets in both the cytoplasm and 
the nucleus 1 A3 For example, in the cytoplasm mRNA 
molecules can undergo what looks like liquid/liquid phase 
separation, to produce droplets enriched in the mRNA 
molecules 17 . But these droplets cannot be at ther¬ 
modynamic equilibrium. They are affected when the 
molecular motors dynein or kinesin are knocked dowrP 3 . 
These motors consume energy and actively move mRNA 
molecule^. This suggests that mRNA molecules do not 
simply diffuse into and out of these liquid droplets, they 
are actively transported into or out of these droplets. 

Inspired by this, we wanted a simple model of density- 
dependent dynamics that generates condensation. These 
dynamics should be via hops from one point to another; 
the hops do not conserve momentum, and there is no 
well defined velocity. This is our simple model of a mo¬ 
tor translating mRNA molecules along microtubules in a 
cell. It turns out that very simple models show conden¬ 
sation. Rather generically, liquids appear whenever the 
hopping rate of a molecule decreases as the local density 
of molecules increases. These out-of-equilibrium dynam¬ 
ics are illustrated in a schematic in Fig. [I] The dynam¬ 
ics stabilise liquid droplets as evaporation of molecules 
from a droplet’s surface is reduced by the slow hopping 
rate of molecules out of the dense liquid. Attractions 
slow hopping rates but any non-equilibrium mechanism 
that also slows hopping out of dense regions will tend 
to have a similar effect, whether or not that are any at¬ 
traction energies directly involved. The Edinburgh group 
and otheraJM 22 l, j lave extensively studied a system with 
rather different microscopic dynamics; in their systems 
the particles have a well-defined velocity. However, de¬ 
spite the differences in the microscopic dynamics, here we 
are studying models in which the mobility decreases as 
density increases, just as they have, and so we see qualita¬ 
tively very similar condensation into liquid-like droplets. 
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FIG. 1. Schematic illustrating out-of-equilibrium dynamics 
that depend on the local density of molecules. Molecules are 
shown as blue discs, the lines with arrows indicate the path 
the top molecule has taken. It slowed down when it was near 
the bottom molecule. 

We start with the familiar 2D lattice gas^HSS. Our 
lattice is L by L lattice sites and has N molecules in 
total. This is a very simple model of mRNA molecules 
in a cell, where we effectively integrate over all other 
molecules, to allow us to explicitly consider only these 
mRNA molecules. The interactions are then of mRNA 
molecules in the presence of all these other molecules, 
and the dynamics are those in the presence of these other 
molecules, including motor proteins. As we integrate over 
these other molecules liquid/liquid separation becomes 
condensation of the one species we explicitly consider. 

We start with the familiar Kawasaki dynamie d 23 * 24 ! 
They of course obey detailed balance. For Kawasaki dy¬ 
namics, time is measured in cycles, and in each cycle L 2 
lattices sites are selected at random, one after another. 
When a site is selected then one of its four neighbour¬ 
ing sites is chosen at random. If one of the pair of sites 
is occupied by a molecule and the other is empty, then 
we attempt to move the molecule from one site to the 
other. This is done as follows. The change in the num¬ 
ber of neighbours of the molecule, An, is computed. If 
An > 0 the move is always accepted as then the energy 
decreases or stays the same, but if An < 0 the move is 
rejected with probability 1 — exp[eAn/fcT]. Here e is an 
interaction energy. 

When the ratio e/kT > i. 7 f| 24 l 25 l, a liquid phase forms. 
The attractions between the molecules cause them to 
condense into a liquid phase. Droplets of liquid nucleate 
and grow until there is a single large droplet coexisting 
with a vapour. 
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FIG. 2. Snapshot of configurations of systems of L = 100 by 
L = 100 sites, with 15% of the sites occupied by molecules. 
This is for our out-of-equilibrium model; od = 2.5 and tm = 
8 , which corresponds to the Ising lattice gas at e/kT = 2.5 . 
We have imposed periodic boundary conditions. 


We now introduce a simple model of dynamics where 
the rate a molecule hops at depends on the local envi¬ 
ronment at the start position of its move. At first sight 
these dynamics may look as though they violate detailed 
balance, but we will show that they do not. The dynam¬ 
ics are as follows. Select a lattice site at random. If it 
is occupied by a molecule, attempt to move it to a ran¬ 
domly selected site within a radius Tm ■ If the site chosen 
is already occupied the attempt is always rejected. If the 
chosen site is empty the move is made with probability 
p = exp[—a£>?z n ], for n n the number of neighbours of the 
molecule in its starting position. Here ap > 0 is a pa¬ 
rameter that couples the dynamics to the local density 
of molecules. In this model the higher the local density 
of molecules the slower the hopping rate of a molecule. 
This could be for a number of reasons, such as neighbour¬ 
ing molecules inhibiting motor transport via mechanisms 
such as direct binding or signalling. Alternatively, even if 
motors are not involved, then a local density-dependent 
slow down could be caused by an out-of-equilibrium pro¬ 
cess such as a phosphorylation/dephosplrorylation cycle 
that modulates binding and so diffusion rates. 

Our density-dependent dynamics map onto the stan¬ 
dard Ising lattice gas, as can be seen as follows. To do 
this we need to consider the transition probability be¬ 
tween state i and state j. p t j ; these states i and j corre¬ 
spond to two positions of a molecule on the lattice. For 
transition between a state i with a molecule surrounded 
by n n (i) neighbours and a state j where the molecule has 
n n (j) neighbours, the ratio of the transition probabilities 
i to j, and j to i are 


Pij_ 

Pji 


exp [-a D n n (i)] 
exp [~a D n n (j) 


exp [~a D [n n (i) - n n (j)]] ( 1 ) 


These satisfy detailed balance, and map to the Ising 
lattice gas with e/kT = an- Thus, for example, our 
model has the usual Ising critical point at = 1.76, 
and a reduced surface tension 7 ' = 7 /kT given by On- 
sager’s expressiorP but with an replacing e/fcT: 7 ' = 
od /2 - In ([1 + exp(—au/ 2 ] / [1 - exp(-a£>/ 2 ]). 

The result of a simulation run at ap = 2.5 is shown in 
Fig-i We see that, as it must, it shows the usual Ising 
lattice gas vapour/liquid coexistence. We fix cm = 8, as 


varying it over a wide range of values (2 to 24) is found 
to have little effect on the dynamics, and cannot change 
the phase behaviour as varying tm does not alter the 
mapping to the Ising lattice gas. 

In terms of studying vapour/liquid coexistence in 
energy-consuming systems, the studies cl osest to this 
work are those of the Edinburgh groupP^ 21 J 22 ) w } 10 
model active Brownian particles with motilities that de¬ 
pend on the local density of these particles, but no attrac¬ 
tions. Redner et a/P and Mognetti et alW *- 1 also study 
systems with local slow down in the dynamics that show 
condensation. In the Edinburgh group’s work, the parti¬ 
cles condense into a liquid phase coexisting with a dilute 
phase. Our models behave in a qualitatively identical 
way, the fact that for the model here high densities de¬ 
cease a hopping rate whereas it decreases a velocity in 
the model studied in Edinburgh appears to make little 
difference. 

However, our motivation is different. The Edinburgh 
group focus on showing that a specific out-of-equilibrium 
model shows equilibrium-like vapour-liquid phase sep¬ 
aration, whereas we are interested in constructing the 
simplest possible models model the liquid droplets in¬ 
side cells. There are also similarities between our model, 
the models studied by the Edinburgh group and others, 
and cooling granular media. In cooling granular media 
locally high rates of inelastic collisions in dense regions 
slow the particles in these dense regions, which in turn 
causes clustering that is partly analogous to condensa¬ 
tion into a liquid-like phase. See Paul and DasP and 
references therein for recent work in this area. However, 
there are differences, these clusters are not at steady state 
and there may well not be a well-defined surface tension 
ther^P 

mRNA-rich droplets in the cytoplasm of cells are 
not the only example of self-organised domains in cells. 
There is a huge amount of self-organised structure in 
the chromatin in the nuclei of cells. Chromatin is the 
DNA of our genes, together with associated proteins and 
RNAlHdlHHUI Examples of this structure are d omains 
called nucleoli where ribosomes are mad^ * 10 ! 11 ! Cajal 
bodies where RNA splicing occurs 2 ^, and what is effec¬ 
tively microphase separation between chromatin t hat is 
being transcribed and chromatin that is not activP 1 - 0 12 . 

So chromatin, which includes huge polymers of DNA, 
exhibits behaviour reminiscent of the microphase separa¬ 
tion seen in block copolymers with immiscible blockPl 
The nucleus is microphase separated in the sense that 
lengths of chromatin that have a common feature, e.g., 
are making ribosomal RNA in the case of nucleoli, have 
separated out from the rest of the DNA. However, the 
dynamics of thejchromatin depends on energy consum¬ 
ing processe d — I— I, unlike in block copolymers where the 
monomers move via thermal diffusion. In conventional 
block copolymers the microphase separation is driven by 
intermolecular attractions between like monomers being 
stronger than between unlike monomers. The mechanism 
may be different for chromatin. 
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FIG. 3. Snapshot of a configuration of a 2D lattice block 
copolymer of length M = 100, with 50 of each monomer type. 
The monomers in green (type D) have dynamics that slow in 
the presence of other monomers, whereas those in red (type E) 
interact only via excluded volume. The parameter ctDP = 0.6. 


Nucl ear s elf-organisation is a very active field of 
researctP2E ii an d we cannot answer all the many un¬ 
solved questions here. However we can see if density- 
dependent dynamics of the monomers of a block copoly¬ 
mer, can result in microphase separation. This would 
be a very simple model of separation between active and 
inactive chromatin, or between a nucleolus and the sur¬ 
rounding chromatin that is not involved in the making of 
ribosomes. 

We introduce a simple out-of-equilibrium block copoly¬ 
mer model. The polymer consists of a linear chain of M 
monomers on a 2D square lattice. The monomers are 
held together by bonds. If the x and y coordinates of 
monomer i are Xi and yi then the bonds along the back¬ 
bone are enforced by insisting that |xj + i — Xi\ < 1 and 
\yi+i — y,\ < 1 for i = 1,M — 1. The monomers are 
of two types: Type E which just interact via excluded 
volume interactions, and type D which in addition have 
dynamics that depend on the local density. 

A simulation cycle consists of M attempted moves, in 
each of which one of the M monomers is selected and an 
attempt is made to move it to one of the eight neighbour¬ 
ing lattice sites. These eight are the sites above, above 
right, right, below right, etc, of the central site. If a type 
E monomer has been selected, that move is always suc¬ 
cessful unless it would move the monomer to an already 
occupied site, in which case it is always rejected. If a type 
D monomer is selected then the move is also rejected if 
it would move the monomer to an already occupied site. 
But if the site is vacant, the move is only made with a 
probability p = exp[— appriMN]- Here umn is the to¬ 
tal number of monomers (excluding the monomers it is 
bonded to) in the eight lattice sites that surround the 
site the molecule would be moved from. The parame¬ 
ter cxdp controls the density dependence of the dynam¬ 
ics. The polymer is simulated on a large lattice; peri¬ 
odic boundary conditions are not used. Note that this 
model for density dependent dynamics also satisfies de¬ 
tailed balance, for the same reason as for our monomeric 
model. The model has hard constraints, both restrict¬ 
ing the molecules to no more than one per site and re¬ 
stricting the bond length between successive monomers. 


For two states i and j that are allowed as they sat¬ 
isfy these requirements the transition probabilities satisfy 
Pij/Pji = exp [-a DP [n M N(i) ~ n M N{j)}}, and so satisfy 
detailed balance. 

For a block copolymer of the two types of monomer, 
the two halves of the polymer separate, with the type D 
monomers forming a condensed globule, while the type 
E monomers form an extended chain. A snapshot of this 
is shown in Fig. [3] 

Note that our model for microphase separation is dif¬ 
ferent to Ganai et aZ.^s model for separation of chro¬ 
matin that is being transcribed from chromatin that is 
not being transcribed. In their work domain forma¬ 
tion is driven by a large increase in effective tempera¬ 
ture in chromatin that is being transcribed, without a 
direct mechanism for density-dependence of the dynam¬ 
ics. Within our model, transcriptionally active regions of 
chromatin separate if the active dynamics of these regions 
are such that their motion slows when in contact with 
other active regions. Similarly, applied to a nucleolus, 
our model predicts that the nucleolus separates from the 
surrounding chromatin due to ribosome-producing chro¬ 
matin slowing its mesoscale motion in the presence of 
other ribosome-producing chromatin — perhaps because 
it is sharing factors or due to direct attractive interac¬ 
tions. 

We have studied very simple models and found that 
molecules whose hopping rate decreases strongly when 
the local density is high, condense into liquid droplets. 
The simplicity of our model, and the fact that condensa¬ 
tion is also seen in very different models in which the ve¬ 
locity decreases when the local density is high 15 16 l 19 t tHl 
suggests that this condensation is quite generic. The 
model is an equilibrium one in the sense that it obeys 
detailed balance, but is for out-of-equilibrium systems 
that consume energy. Thus to the extent that it mod¬ 
els the microscopic dynamics of mRNA molecules in cells 
correctly, their behaviour will be qualitatively that of an 
equilibrium liquid. Although the model maps to the Ising 
lattice gas, an equilibrium model, the energy driving the 
dynamics does not have to be kT. For mRNA molecules 
in cells, a better candidate, at least for long distance 
motiorfS', is the much larger forces that molecular mo¬ 
tors can provide. This has consequences. For example 
the effective surface tension that limits the fluctuations 
of the liquid surface is then not ~ kT/a 2 , it is presum¬ 
ably ~ wm/o : 2 , with u>m the work done by one or a few 
motors on one of the molecules of the liquid. Here a is a 
molecular diameter. 
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